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Abstract. Gamma-ray bursts are believed to be produced in highly-relativistic collimated outflows. Support for 
this comes among others from the association of the times of detected breaks in the decay of afterglow light curves 
with the collimation angle of the jets. An alternative approach to estimate a limit on the collimation angle uses 
GRB afterglows without detected prompt-emission counterparts. Here we report on the analysis of a dedicated 
survey for the search of these orphan afterglows using the Wide Field Imager at the 2.2 m MPI/ESO telescope at 
La Silla, Chile. We monitored ~12deg 2 in up to 25 nights typically spaced by one to two nights with a limiting 
magnitude of _R=23. Four previously unknown optical transients were discovered and three of these associated 
with a flare star, a cataclysmic variable and a dwarf nova. The fourth source shows indications for an extragalactic 
origin but the sparse sampling of the light curve prevents a reliable classification. We discuss the results in the 
context of the collimation of GRBs. 
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1. Introduction 

There is now conspicuous observational and theoretical 
evidence that the radiation of cosmic gamma-ray bursts 
(GRBs) is produced in highly relativistic collimated out- 
flows. A jet geometry was originally invoked as a solution 
for the "energy crisis" by reducing the total energy output 
of a GRB by a factor of f2 7 /47r, where J7 7 is the solid angle 
into which gamma-rays are emitted (e.g., Rhoads 1997; 
Fruchter etal. 1999). Evidence for the collimation was 
provided by the theoretically predicted (e.g., Sari 1999; 
Gruzinow 1999; Ghisellini & Lazzati 1999) and observed 
polarization evolution of optical afterglows (e.g., Covino 
etal. 1999; Wijers etal. 1999; Greiner etal. 2003). 

Prominent signatures of the jet geometry were identi- 
fied already earlier in the broad-band breaks observed in 
the optical and radio afterglow light curves of several long 
duration GRBs (e.g., Stanek etal. 1999; Harrison etal. 
1999). The opening angles, which have been inferred from 
these "jet breaks" for a number of afterglows, vary from 
1° to more than 25°, with a strong concentration near 4° 
(Frail etal. 2001). 

A consequence of the collimation is that the prompt 
7-ray emission will be detected if the viewing angle of 
the observer is equal to or smaller than the opening an- 
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gle of the jet, 6j e t- This implies that only for a fraction 
of all bursts in the Universe 7-ray photons will reach the 
Earth. The total GRB rate will be higher than the ob- 
served rate by a factor of roughly 9~ 2 t oc V 2 in case of a 
universal jet structure, where V is the bulk Lorentz factor 
of the ejecta. For a quasi-universal Gaussian-type jet, the 
rate will be generally smaller. Nevertheless, GRBs which 
are not pointed directly at the observer can in principle 
be discovered through their afterglow radiation at longer 
wavelengths. 

In the standard internal-external fireball model (e.g. 
Rees & Meszaros 1992, 1994) the afterglow is produced 
when the initially highly relativistic ejecta plows into the 
ambient medium. The afterglow emission is radiated into a 
solid angle of fi^~l/r, along the line of motion. When the 
jet decelerates, fi^ increases until it reaches 4tt. Therefore, 
so-called "off-axis" orphan afterglows can be detected for 
bursts which are beamed outside of the field of view of the 
observer (e.g., Rhoads 1997; Perna & Loeb 1998; Dalai 
etal. 2002). The light curves of these off-axis orphan af- 
terglows will initally be faint, brightening up to a viewing 
angle dependent maximum and become similar to regular 
GRB afterglows after the jet break later on (Rhoads 1999; 
Nakar & Piran 2003). 

It is theoretically reasonable to assume that the optical 
afterglow might be emitted from a slower moving mate- 



2 



A. Rau et al.: Constraining the GRB Collimation 



rial (a lower T) than the initial 7-ray beam. Therefore, 
so-called "on-axis" orphan afterglows are expected when 
the narrow 7-ray emission misses the observer by a small 
amount but the wider optical emission region falls within 
the observation cone (Nakar & Piran 2003). These after- 
glows will exhibit similar light curves as regular afterglows 
with detected prompt emission. 

Observations of orphan afterglows can help to study 
the initial opening angle of the jets and to place a con- 
straint on the collimation of the optical afterglow emission 
(Rhoads 1997). Especially on-axis orphans are suitable as 
they are substantially brighter than off-axis orphans and 
thus easier to detect in a dedicated survey (Nakar & Piran 

2003) . Additionally, Dalai et al. (2002) pointed out, assum- 
ing an uniform jet with constant jet-break time and lumi- 
nosity at the break time for an on-axis observer, that for 
small angles the afterglow effective beaming angle scales 
with the jet opening angle. Therefore, the number of de- 
tectable off-axis orphan afterglows would be independent 
of the jet opening angle and similar for moderately wide 
jets (~20°) and for arbitrarily narrow jets (<0.01°). In 
contrast, Totani & Panaitescu (2002) predict a strong de- 
pendence of the orphan rate on the jet opening angle as- 
suming a constant total energy in the afterglow jet. 

A small number of surveys dedicated to the search 
of untriggered optical GRB counterparts were performed 
over the past years. No candidate event was found in 
125 hrs monitoring of a field of 256 deg 2 with ROTSE- 
I to a limiting magnitude of 15.7 (Kehoe etal. 2002). 
Vanden Berk et al. (2002) searched for color-selected tran- 
sients within 1500 deg 2 of the Sloan Digital Sky Survey 
(SDSS) down to R=19 and found only one unusual tran- 
sient which was later identified as a radio-loud AGN ex- 
hibiting strong variability (Gal- Yam etal. 2002). The au- 
tomated RAPTOR wide-field sky monitoring system al- 
lows to image 1300 deg 2 at a time down to a magnitude of 
^12.5 (Vestrand etal. 2004). A couple of interesting opti- 
cal transients were found in the B, V and i?-band Deep 
Lens Survey transient search, within an area of 0.01 deg 2 
yr with a limiting magnitude of 24. None of those could be 
positively associated with a GRB afterglow (Becker et al. 

2004) . Recently, Rykoff etal. (2005) performed a search 
using the ROTSE-III telescope array without detecting 
any candidate afterglow events. They placed an upper 
limit on the rate of fading optical transients with quies- 
cent counterparts dimmer than ~20th magnitude of less 
than 1.9 deg -2 yr -1 . 

In this paper we present the results of a survey for op- 
tical orphan afterglows with a wide-field imaging instru- 
ment and will give an estimate of the GRB collimation 
derived from the detection rate of on-axis orphan after- 
glows. The paper is structured as follows. In Sect. 2 we 
describe the observational strategy, instrumentation and 
selected survey fields as well as the data reduction and 
transient detection method. The candidate transients are 
presented in Sect. 3 and the transient detection efficiency 



in Sect. 4. A discussion of the results in the context of 
GRB collimation is given in Sect. 5. 

2. Observations and Data Reduction 

2.1. Strategy, Instrumentation &. Survey Fields 

A search for GRB afterglows or other transient phenom- 
ena in the sky requires a thorough strategy due to the 
random occurence of the events in space and time. As a 
continuous monitoring of a large field with a large aper- 
ture telescope could only be considered overambitious in 
many cases, we developed our strategy using the knowl- 
edge of the properties of observed optical afterglows at the 
time when the survey was proposed (1999). 

The primary idea was to take multiple deep observa- 
tions of a number of selected sky fields. The observing 
scheme was chosen such that the survey would be sen- 
sitive enough (i?^23mag) to provide the detection of a 
GRB orphan afterglow in at least two epochs together 
with earlier and later upper limits. To combine the avail- 
ability of a large aperture telescope together with the ob- 
served brightness decay of GRB afterglows, we decided to 
perform consecutive observations of a given field in every 
2nd night. Over this time span, on-axis orphan afterglows 
will in most cases be brighter than R=2\. Instead, off- 
axis orphans are expected to be fainter and will start to 
dominate below that magnitude (Nakar & Piran 2003). 

We obtained imaging data during three periods (May- 
June, August & September-October) in 1999 and moni- 
tored 7 different sky fields in up to 25 nights each. Each 
field is composed of 2 to 8 separate sub-fields and a to- 
tal of 38 sub-fields were selected. This corresponds to an 
area of ~12 square degrees being imaged throughout the 
campaign. A list of all fields together with the number of 
sub-fields and the maximum number of observing nights 
is provided in Tabled Note that due to changes of the ob- 
serving conditions during a night not always a complete 
monitoring of a given field could be accomplished. 

A total of 39 nights were scheduled using the Wide 
Field Imager (WFI) at the MPI/ESO 2.2 m telescope 
in La Silla, Chile. Due to weather constraints only 31 
nights could at least partly be used for observations. 
Unfortunately, the lost nights caused unfavorable inter- 
ruptions of the otherwise gap-less observing schedule of 
each period and lead to a significantly decreased detec- 
tion sensitivity for orphan afterglows. 

The distribution of the time delay between two consec- 
utive observations of a given sub-field is shown in Figure^ 
While the majority of the survey could be observed with 
the proposed gap of two days between consecutive obser- 
vations, ~15% suffered from larger gaps due to inappro- 
priate weather conditions. In addition to the regular 2-day 
schedule a small number of nights were included in which 
multiple observations of selected fields (Fl, F3 & F5) were 
performed within one night. This would allow to identify 
and distinguish possible short-term variable sources (e.g. 
CVs) in these fields. 
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At [days] 

Fig. 1. The distribution of time between consecutive ob- 
servations of individual sub-fields with 1 hr binning. 

The photometry was taken with WFI mainly in the 
i?-band together with a small number of pointings per- 
formed using the V and J-band filters. The typical expo- 
sure time was 420 s per pointing. The instrument consists 
of a mosaic of 4x2 CCDs each 2046x4128 pixel in size 
which, together with the plate scale of 0.238 arcsec/pixel, 
provided a sky coverage of 34' x 33' per image. A gap 
of ~2 arcsec exists between neighboring CCDs of the im- 
ager. Throughout the data acquisition the observing con- 
ditions varied strongly (see Figure^. The majority of the 
data were taken at a seeing of ~1" but a small fraction 
of the observations suffered from significantly worse see- 
ing. Due to the low limiting magnitude and the loss of all 
but the brightest sources we neglected date taken under 
seeing conditions above 2'.'4 (~10 pixel FWHM). This cor- 
responds to an additional loss of ~10% of the campaign 
data. 



0.5 1 1.5 2 2.5 

seeing [arcsec] 

Fig. 2. Seeing distribution for all pointings. 

The selection of the fields was more or less random 
with respect to the search for orphan afterglows with the 
exception of the field F6. This field is centered on Selected 
Areas 113 and includes a number of Landoldt standard 
stars which allow an absolute photometric calibration of 
the survey (Landoldt 1992; see below). F6, together with 



Table 1. Observation log. The first four columns provide 
the ID, central coordinates and number of sub-fields. The 
last column represents the amount of nights with seeing 
<2'.'4 in which pointings of a given field were obtained. 



Field 


RA(2000) 


Dec(2000) 


sub-fields 


# of nights 


Fl 


01h32m 


-43° 12' 


4 


15 


F2 


03h33m 


-27° 37' 


4 


12 


F3 


13h28m 


-21° 40' 


8 


11 


F4 


16h20m 


+04° 00' 


8 


12 


F5 


21h26m 


-43° 22' 


8 


23 


F6 


21h41m 


+00°30' 


2 


25 


F7 


21h52m 


-27°32' 


4 


21 



a second field (F4), is also covered by the SDSS 1 . This 
provided us with the opportunity to identify possible tran- 
sient sources in these fields down to i?~22mag using the 
SDSS multi-band informations. Particular care was taken 
to avoid stars brighter then ~12th magnitude which would 
not only saturate the detector but also leave effected pix- 
els less sensitive for subsequent images. Furthermore, high 
ecliptic lattitudes were favored to maximize the Moon dis- 
tance, and thus make the fields observable over a large 
fraction of an observing period. Finally, significant fore- 
ground extinction was avoided during the field selection. 



2.2. Data Reduction 

Nearly 700 images were obtained throughout the survey 
comprising ~130GB of raw data. For the automatic re- 
duction and analysis of this large amount of data a Perl- 
based pipeline was developed which uses a number of well 
tested astronomical software packages. The basic image re- 
duction was performed using IRAF 2 /MSCRED. For each 
night of observations a common bias frame was produced 
and subtracted from the science images. Flat field correc- 
tion was performed using daily super-sky-flats produced 
from all science observations of a given night without sig- 
nificant illumination by the Moon or other bright sources. 

The astrometric solutions were obtained separately for 
each of the 8 CCDs with the WIFIX/ ASTROMETRIC 
package by comparing the positions of detected sources 
with those compiled in the USNO-A2.0 catalog (Monet 
etal. 1998). The resulting astrometric precision is indi- 
cated in Figure |3 where the difference in the position of 
isolated sources detected in a sub-field of F6 in multi- 
ple observations is shown. About 80 % of the sources lie 
within 1/2 detector pixel (0'.'119) compared to a refer- 
ence observation taken on May 31, 1999. More than 99% 
of the sources are detected within a circle of 0'.'4 radius. 
Note that the double peaked shape of the distribution in 



1 |http:/ / www.sdss.org/ 1 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 

3 http:/ /www. na.astro.it/~radovich/wifix. htm 
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Figure|3]is artificial, resulting from the numerical rounding 
of the pixel coordinates in the source detection algorithm. 



0.2 



0.2 0.4 0.6 0.8 

position difference ["] 

Fig. 3. Normalized distribution of the distances of ~40 
isolated sources detected during the whole survey in the 
field F6-1 compared to the respective reference positions 
obtained from an observation taken on May 31, 1999 (solid 
line). The dotted and dashed lines indicate the size of 1/2 
detector pixel and the search radius for the source tracing 
algorithm, respectively. The selected sources span ^5 mag 
in brightness and are distributed over the entire detector 
field of view. 

The source detection in the images was performed 
in each of the 8 WFI CCDs separately. We used the 
IRAF/DAOPHOT package to measure the source flux in- 
side a Gaussian-shaped point-spread-function. The search 
for transient sources does not gain significant benefit from 
an absolut photometric calibration which would be con- 
nected with an unavoidable additional systematic zero- 
point uncertainty. Instead, we obtained the variability 
information using the technique of differential photome- 
try. We selected an ensemble of at least 20 local, non- 
saturated, non-variable reference stars for each sub-field. 
By deriving the median brightness offset of these stars 
with respect to the brightness of the same stars obtained 
in a reference image of a given sub-field (typically the one 
with the smallest seeing), the photometric offset between 
the two observations could be estimated. This was suc- 
cessively done for all observations of a given sub-field and 
thus provided a common photometric zero-point for all 
pointings of the respective sub-field. 

Absolut photometric calibration was obtained based 
on observations in one photometric night (October 12, 
1999) of the field SA113 (our F6). This field contains a 
number of Landoldt standard stars with well tabulated 
optical photometry (Landoldt 1992) 4 . As all photometric 
standard stars were saturated in the regular 420 s image 
taken for our survey, we obtained a shorter exposure (20 s) 



during the same night. In this 20 s frame six non-saturated 
standard stars were contained and the photometric zero 
point for this image was obtained (systematic uncertainty 
of Ai?=0.1 mag). This allowed to produce a sample of sec- 
ondary standard stars in the field down to a limiting mag- 
nitude of i?^21 corresponding to a 5a detection. 

The field F6 was observed in nearly all nights of the 
survey. This allowed to obtain accurate photometric cal- 
ibration for the other six fields based on the zero points 
derived from the secondary standards in F6 in a common 
night. For each field one photometric night was chosen 
and the photometric calibration of all further observations 
were calibrated with respect to this night. The photometry 
in all fields was additionally corrected for Galactic fore- 
ground extinction (,4 fl =0.02-0.25) (Schlegel etal. 1998). 

Figure 0] shows the resulting photometric quality for 
an example field (F4) obtained under good seeing condi- 
tions (078). The uncertainties represent the detection un- 
certainties in the instrumental magnitude system and do 
not include the uncertainties resulting from the absolute 
photometric calibration (Ai?~0.1 mag). A limiting mag- 
nitude of i?~23 at 10a was reached in this observation. As 
discussed above, the observing conditions varied strongly 
throughout the survey. Thus, the achieved limiting mag- 
nitude ranged between i?^19 and i?~23, depending on 
the seeing and Moon illumination. Sources brighter than 
i?~16 were saturated in a good night. 



0.15 




4 The tabulated Cousin 7?-band magnitudes were trans- 
formed to the WFI 7?-band filter system using the colour and 
extinction terms given at . 



J i i i I i i i L_ 

18 20 22 

R-band Brightness [mag] 



Fig. 4. Photometry obtained for an observation of the 
field F4 on May 22, 1999 with each dot representing one 
source. The seeing during the time of the imaging was 
Of! 8 and a limiting magnitude of i?~23mag with lOcr was 
reached. Sources brighter than i?^16mag are saturated. 



2.3. Transient Detection Method 

For each of the seven fields a master table including the 
coordinates and magnitudes of all detected sources in all 
observations was produced. Herein, we assigned detections 
within a distance of 0'.'75 to a detection obtained in a 
reference image together and consider this as one source. 
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Light curves were obtained for all detected sources and 
candidate transient objects were selected based on the de- 
viation of their light curve from their mean light curve. 
For all candidates with single detections or Ai?>0.75mag 
(^12000), light curve plots as well as thumbnail images 
of all pointings were produced. These were examined by 
eye in order to remove spurious transients arising from 
nearby bright stars, extended objects, stray light effects 
at the edge of the FoV, bad focus or detections which 
were consistent with faint stars at the limiting magnitude 
of the individual pointings. The strategy of the survey was 
aimed to catch a possible orphan afterglow in at least two 
consecutive observations. Therefore, we considered only 
sources with at minimum two detections for further inves- 
tigations. A number of fast moving objects (e.g. air planes, 
terrestrial satellites) could be identified by the trace in the 
images which they left during the 420 s exposure. Solar 
system objects (e.g. asteroids) are especially abundant in 
the survey fields located close to the ecliptic plane (F3, 
F6 & F7). These sources appear either as one-time de- 
tections or show significant motion between subsequent 
observations and could be identified accordingly. Of the 
~12000 candidates only four remained and were double 
checked with the positions of known sources compiled in 
the SIMBAD 6 , NED 7 and NSV 8 databases and if possible 
correlated with the SDSS 4th release and the DSS 9 . 



3. Results 

3.1. Candidate Transients 

Throughout the survey four new transient sources with 
detections in at least two images were found. Based on 
the shape of the light curve and amplitude we identified 
one candidate cataclysmic variable, one dwarf nova, one 
flare star and one candidate extragalactic transient su- 
perimposed on an underlying faint object. For the latter 
an orphan afterglow nature is considered plausible. Below 
we give a brief description of each of these four transient 
sources (see also Table EJ. 



Table 2. Candidate transient sources. The minimum and 
maximum obtained brightnesses as well as the putative 
classification are given. 



Name 



Brightness [mag] 
min max 



Putative 
classification 



J132653.8-212702 
J132813.7-214237 
J161953.3+031909 
J215406.6-274226 



20.3±0.1 
21.3±0.1 
19.9±0.1 
22.5±0.2 



19.5±0.1 
19.9±0.1 
17.5±0.1 
20.0±0.1 



CV, eclip. binary 
extragalactic 
dwarf nova 
flare star 



6 http://simbad.u-strasbg.fr/ 

7 http:/ /nedwww. ipac.caltech.edu/ 

8 http:/ /heasarc. gsfc.nasa.gov/W3Browse/all/gcvsnsvars. html 

9 http:/ /archive. stsci.edu/cgi-bin/dss_form 



3.1.1. J132653.8-212702 

The source was initially detected in F3 during the first 
night of observations on May 22, 1999 (MJD 51321.0045). 
The field was observed four times during that night and 
the source exhibited a significant brightness increase from 
i?=20.4 to i?=19.7 within 80min (Fig. Et). ft was ob- 
served again on May 31/ June 2, June 17/19 and August 
4/6 and showed a flaring of Ai?=0.2-0.4mag during each 
of these epochs. No X-ray or optical counterpart was de- 
tected in the ROSAT All-Sky survey (taken 1990). No 
entry was found in SIMBAD or NED. The rapid variabil- 
ity excludes slowly varying objects and is consistent with 
a cataclysmic variable and/or eclipsing binary. No orbital 
period could be obtained due to the sparse sampling of the 
light curve. The confirmation of the nature of this source 
requires further monitoring of the variability or simultae- 
nous observations in multiple colours. 



3.1.2. J132813.7-214237 

The source was detected on June 20.07 UT 1999 
(MJD 51349.0430) at its maximum brightness of J?=19.9 
(Fig. EtO- Unfortunately, the observing period finished 
after the observation of the outburst and the next i?- 
band pointing was performed more than six weeks later 
on August 04. Thus, a decay of the source could not be 
monitored. A faint (i?=21.3) persistent point source was 
found at the position of the transient in images taken 
before and after the outburst. The ROSAT All-sky sur- 
vey did not show a source at the position of J132813.7- 
214237 at 3cr flux upper limit of 5.5x 10~ 13 erg cm' 2 s -1 . 
Similarly, no counterpart was found in the SIMBAD and 
NED databases. The object is visible near the limiting 
magnitude of the DSS in the B and i?-band and not de- 
tected in the DSS /-band. A detailed analysis revealed 
that the transient source was offset by ~0>.'8 from the po- 
sition of the persistent counterpart (Fig. EJ). This might 
suggest an extragalactic origin of the transient assuming 
the persistent counterpart can be associated with a can- 
didate host galaxy. Nevertheless, this vague assumption 
based on the available observational data requires a con- 
firmation by an accurate distance measurement. 

The lacking observational coverage of the decay light 
curve of J132813. 7-214237 leaves the identification of the 
=origin of the source an open question. Assuming the ex- 
traggalactic nature of the faint persistent source, the flar- 
ing source could, for instance, be associated with a su- 
pernova explosion occuring in this galaxy. The bright- 
ening of ~1.5mag in ^2 days is very steep compared 
to observed supernovae though (Leibundgut etal. 1991). 
"Another explanation is that of a foreground flare star close 
to the line of sight towards the persistent background 
source. Furthermore, J132813. 7-214237 might be a poten- 
tial orphan afterglow. In order to test this hypothesis we 
searched for triggered GRBs which occurred during the 
time between the preceding observation (June 18.05 UT) 



6 



A. Rau et al.: Constraining the GRB Collimation 




51320 



51400 



(b) 



t 

H — h 



51320 



51340 



51360 
MJD 



H — I- 



513B0 



51400 




51350 



51400 



_ 1 1 1 1 








■ 












■ 












I 






20 


z 1 

r ■ 


hi" 








20.2 












20.4 




■ 








20.6 






















- + 


20.8 

















51343.42 


51343.43 




: i i i 1 


1 1 1 1 h 






■| — 1 — 1 — 1 — I 
■ 












■ 




— ■ 












■ 

— T T 1 1 













51330 



MJD 



51335 
MJD 



51340 



Fig. 5. Light curves of the four transient source candidates (top) and seeing evolution (bottom). Error bars are 
smaller than the symbol size if not visible, top left: The candidate cataclysmic variable J132653. 8-212702. top right: 
J132813. 7-214237, extragalactic transient candidate, bottom left: J161953. 3+031909, a dwarf nova, bottom right: 
J215406. 6-274226, a possible flare star or afterglow candidate. The inset shows the detailed light curve at the time of 
the flare. 



and the outburst. We found 5 cataloged 10 GRBs during 
this period (BATSE #7609 & 7610 and IPN #2066, 2067 
& 2069; K. Hurley, private communication). None of those 
has a position consistent with J132813. 7-214237 which 
allows to exclude an association with a triggered GRB. 
Regardless, J132813. 7-214237 can also be due to an un- 
triggered or orphan afterglow and is the best such candi- 
date found during our survey. 

3.1.3. J161953.3+031909 

This candidate was detected as a constant source 
(i?=19.9) during the first four observations of the field F4 
in May and June 1999 and exhibited a sudden brighten- 
ing by Ai?=2.4mag between the pointings on June 14.29 
UT (MJD 51343.2950 ) and June 17.14 UT (Fig. 0:). 
The following observation on June 19.15 UT showed the 
source unchanged and later observations indicate a subse- 

10 http:/ /grbcat. gsfc.nasa.gov/grbcat/grbcat. html 



quent decay over 50-90 days back to the quiescent bright- 
ness. The optical light curve suggests the classification of 
J161953. 3+031909 as a dwarf nova. This is strengthened 
by the detection of a faint X-ray source in the ROSAT 
all-sky survey. During an exposure of 320 sec on August 
12/13, 1990 a total of 9 source photons were detected, 
corresponding to a mean vignetting-corrected count rate 
of 0.033 cts s _1 (this is below the significance threshold 
of the all-sky survey catalog, so this source is not con- 
tained in the 1RXS catalog of Voges et al. 1999). Adopting 
a thermal bremsstrahlung model with 1 keV temperature 
and half of the Galactic foreground absorption, an unab- 
solved flux in the 0.1-2.4keV band of (7.0+1.0) x 10~ 13 erg 
cm~ 2 s _1 (or (9±2)xl0~ 13 erg cm -2 s _1 bolometric) is 
derived. Using the quiescent optical brightness, this im- 
plies a ratio of Lx/L opt = 0.6, consistent with SU UMa 
stars (Verbunt et al 1997). With an X-ray luminosity of 
1.1 xlO 32 [D/l kpc] erg s^ 1 , the implied distance is of order 
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Fig. 6. Finding charts for the four discovered transients. From top left to bottom right: J132653. 8-212702, J132813.7- 
214237, J161953.3+031909, J215406.6-274226. North is up and East to the left. 



a few hundred parsec. No entry was found in the SIMBAD 
and NED databases. 

3.1.4. J215406.6-274226 

This source was initially detected in two epochs of the field 
F7 as a faint object with i?=22.5 (Fig.EJl). On June 14.41 
UT 1999 (MJD 51343.4150) it was found Ai?=1.9mag 
brighter than previously measured. Fortunately, during 
this night four images of the field were taken and the rapid 
fading of J215406.6-274226 by Ai?=0.7mag in around 
20 min, corresponding to a decay with t , was discovered 
(inset of Fig. . Unfortunately, the further fading could 
not be monitored as the position of the source fell into the 
gap between two CCDs for the following twelve pointings. 
While the rapid decay would be consistent with the obser- 



vations of early GRB afterglows, also a flare star offers a 
possible explanation. Given the range of absolute magni- 
tudes for nearby M-dwarfs of M v =12-16 (Reid et al. 1995) 
the observed quiescent brightness would place J215406.6- 
274226 at a distance of 0.2-1.2 kpc. No X-ray or optical 
counterpart was detected in the ROSAT All-Sky survey 
and SIMBAD or NED, respectively. No source is detected 
in the DSS B and i?-band images but a faint source (^7cr) 
is visible in the DSS /-band. The quiescent counterpart 
could not be resolved and further observational effort is 
required for a more solid classification of the object. 

4. Efficiency for Orphan Afterglow Detection 

The detection efficiency for on-axis optical afterglows was 
estimated using a set of Monte-Carlo simulations folded 
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Fig. 7. left: Image of J132813. 7-214237 (center source) 
and a nearby bright star from an observation taken on 
May 23.00 UT when the source had a brightness of R=2\. 
Flux contours are overplotted for both objects. North is 
up and East to the left, right: Same field observed on 
June 20.07 UT together with the flux contours from May 
23.00 UT. The flaring source shows an offset of ~0'/8 with 
respect to the quiescent counterpart. 



with the observing schedule of the survey. We simulated 
afterglows with random sky coordinates, light curve pa- 
rameters and explosion times distributed over the periods 
in which observations were taken. For all afterglows with 
positions inside one of the monitored fields, the expected 
magnitudes of the afterglows in the first and second ob- 
servation of the field after the burst were calculated. 

The afterglow light curves were described by a broken 
power law and parametrized by the pre-break slope, a±, 
break time, tf,, post-break slope, and initial i?-band 
magnitude, Ri n . Here, the initial magnitude corresponds 
to the observer frame i?-band brightness with which an af- 
terglow is created in the simulations. We used the param- 
eter ranges of 38 observed optical afterglow light curves 
summarized in Zeh et al. (2005). A relatively flat decay 
plus a steeper flux decrease after the jet break is typically 
observed for optical afterglows associated with detected 
long-duration GRBs and is also expected for on-axis or- 
phans as well as "regular" afterglows of untriggered bursts. 
As described earlier, off-axis orphans will show a different 
behavior. They are expected to be fainter than on-axis or- 
phans at early times and follow the post-break decay of on- 
axis afterglows after an initial phase of re-brightening. It 
appears safe to assume that on-axis afterglows will be the 
majority at a limiting magnitude of i?=19-21 and off-axis 
orphans will dominate only at lower magnitudes (Nakar 
& Piran 2003). As the survey strategy foresees to consider 
only sources as candidate orphans with at least two de- 
tections spaced by two nights, off-axis orphans with max- 
imum brightness of i?>21 will in most cases be to faint to 
be identified. While bright off-axis orphans could in prin- 
ciple have been detected in the survey, we do not include 
them in the simulations and focus solely on untriggered 
and on-axis orphans. 



The expected number of detected afterglows in the 
survey depends on the choices for a\, a.2, tb and Ri n - 
The ranges for the slopes and break time were obtained 
from observed light curves (Zeh etal. 2005). Accordingly, 
we used Gaussian distributions for -1.8<ai<-0.4 and - 
2.8<a2<-1.4 and an uniform distribution for tb ranging 
from 0.4-4 days. The most influencing parameter for the 
outcome of the simulations is the initial magnitude of 
a candidate afterglow. Rapid optical follow-up observa- 
tions of BATSE bursts with LOTIS and ROTSE-I (Park 
etal. 1999; Akerlof etal. 2000; Kehoe etal. 2001) and re- 
cent Swift/XJVOT detections indicated that the prepon- 
derance of early afterglows does not get brighter than 
R~U. Furthermore, HETE-2 follow-up and UVOT obser- 
vations showed that around 50 % of the afterglows might 
be brighter than i?~18.5 for 30 minutes after the burst 
(e.g., Lamb etal. 2004). 

In order to test the influence of the limits of Ri n on the 
expected number of afterglows in the survey, we performed 
two simulations with 10 6 bursts per year and full sky, 
each. The initial magnitudes were uniformly distributed 
between 9<i?i„<20 and 13<i?i„<23, respectively. In addi- 
tion, a distribution proportional to 0.2xi?i„ in the range of 
13<i?i ra <23 was simulated. The latter corresponds to the 
observations of afterglows presented in Zeh et al. (2005) 
and is expected to reproduce the reality more closely than 
an uniform distribution. Each simulation was repeated 10 3 
times and the mean detection rates at a given i?-band 
magnitude were obtained. 

The results of the simulations are presented in 
Figure |S] Here we show the normalized (devided by the 
total number of simulated events) probability for a first- 
time detection of a simulated afterglow brighter than a 
given magnitude. The normalized probability represents 
the chance to detect the afterglow in a simulation with 
one single event over the year and full sky. The expecta- 
tion for a second observation of a candidate afterglow for 
the model with non-uniform distribution of Ri n is included 
as well. Naturally, the probability to detect an afterglow 
increases with the depth of the survey. More transients 
are expected to be detected above a certain magnitude 
for afterglow models with brighter R in . Due to the lower 
number of bright afterglows and the increasing probability 
with diminishing initial brightness, the non-uniform dis- 
tribution gives the lowest expectation rates. At the limit- 
ing magnitude of the survey of i?~23 the probability to 
detect a specific randomly selected afterglow in at least 
one observation is approximately 3xl0 -7 . The detection 
rate of events in at least two consecutive observations is 
lower by a factor of ~3. In the more optimistic models, the 
probability for a single detection reaches up to 2xl0~ 6 . 

5. Discussion 

The Monte Carlo simulations described in the previous 
section provide us with the number of afterglows (per year 
and full sky), Nmc, which correspond to the probabil- 
ity of identifying one event in two consecutive observa- 
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Fig. 8. Probability distributions for the detection of a ran- 
domly selected afterglow for three different models of the 
initial brightness. The dotted and dashed functions repre- 
sent uniform distributions of i?i„=9-20 and i?i„=13-23, 
respectively. The thick and thin solid lines correspond to 
the first-time and second-time observations for a distribu- 
tion proportional to 0.2xi?i„ in the range of i?i„=13-23. 
The vertical dotted line marks the limiting magnitude of 
the presented survey. 

tions of the survey. For the three tested afterglow param- 
eterizations we find values between iVjv/c=l-5xl0 6 and 
iVA/c=lxl0 7 . 

Throughout the survey, four unidentified optical tran- 
sients were discovered and one of these sources shows 
indications for an extragalactic origin (see Sect. 3.1.2). 
Nevertheless, the flaring of the source was only detected in 
a single observation and thus the fading could not be moni- 
tored. An unambiguous identification of the transient with 
a GRB afterglow is therefore not possible which leaves us 
with the result of having no clear orphan or untriggered 
afterglow detected in the data. Therefore, Nmc obtained 
from the simulations can be interpreted as an upper limit 
on the true number of on-axis afterglows per year and full 
sky. For N>Nmc one or more on-axis orphans would have 
been expected in the data. 

As advertised earlier, on-axis orphan afterglows can be 
used to place a constraint on the collimation of the opti- 
cal afterglow emitting region relative to the collimation 
of the 7-ray emitting jet. The collimation factor, / c , cor- 
responds to the ratio of the true rates of on-axis optical 
afterglows, Na, and long-duration GRBs which produce 
observable optical afterglows, iV 7 , pointed at the Earth. 
With Nmc>Na, an upper limit for the collimation can 
be derived from f c <NMc/N 1 , Obtaining iV 7 is not triv- 
ial as several uncertain factors influence the number of 
GRBs for which afterglows are in principle detectable. In 
a simplified model iV 7 can be written as 

iV 7 = N 7 , obs -fx-fo-fs (1) 

where N ly0 b s is the observed rate of long-duration 
GRBs with a specific 7-ray instrument corrected for sky 



coverage, Earth-blockage and instrument down-time fx 
corrects for events outside the instrument's energy range, 
fo is the correction factor for optically dim or dark af- 
terglows of long-duration bursts (e.g., intrinsically faint, 
absorbed, high-z) and fs corrects for possible afterglows 
associated with short bursts. Using the full sky GRB 
rate measured with BATSE of ~666 yr _1 (Paciesas etal. 
1999) and correcting for the ratio of long to short bursts 
(2:1; Kouveliotou etal. 1993), we obtain A^ 7 , obs — 444. The 
HETE-2 population of high energy bursts revealed a com- 
position of X-ray flashes (Heise etal. 2001), X-ray rich 
bursts and "normal" GRBs in equal parts (Lamb et al. 
2005). X-ray flashes and X-ray rich bursts show similar 
afterglows as observed for GRBs but have softer prompt 
emission spectra. As BATSE was less sensitive at lower 
energies than HETE-2 a fraction of these events might 
have been missed. This could require a correction factor 
as large as /x=2. 

In general, not for all rapidly followed GRBs an as- 
sociated optical transient can be found. Some detected 
afterglows were faint already early on and would have 
fallen below the limiting magnitude reached in our sur- 
vey. The fraction of these events was found to be of the 
order of 10 % (e.g., Lamb et al. 2004; Jakobsson et al. 2004; 
Rol etal. 2005) which provides /d^0.9. Only recently the 
first optical transients of short-duration GRBs were found 
(Price et al. 2005; Gal- Yam et al. 2005). Theory and obser- 
vations hint that these afterglows are significantly fainter 
than the counterparts of long-duration bursts. Therefore, 
we neglect the influence of short burst afterglows and ap- 
ply /s=l (as we also subtracted them from the burst rate 
of BATSE). 

Using the assumptions discussed above together with 
AT MC =lxl0 7 we derive iV 7 ~800 and / c <12500 accord- 
ingly. This rather conservative upper limit is significantly 
higher than the beaming correction derived by Guetta 
etal. (2005) and Frail etal. (2001) of 75±25 and 500, re- 
spectively. The 7-ray beaming factor corresponds to the 
ratio of the overall rate of GRBs to the detected burst 
rate and should be an upper limit for fc (corresponding 
to the case of isotropic afterglow radiation). 

The high upper limit on fc shows that the effective 
coverage of the performed observations were not sufficient 
to provide a strong constraint for the collimation. An ap- 
proximately 25 times larger portion of the sky (~325 deg 2 ) 
would have been required to reduce f c to less than 500, 
assuming no orphan afterglow detection. We performed a 
further Monte Carlo simulation in order to estimate the 
properties for an "ideal" survey assuming a schedule of one 
observation per field every two nights over 150 nights with 
a limiting magnitude of i?=23. We find that /c<500 (<75, 
<10) would be reached with such a configuration and a 
50 deg 2 (300 deg 2 , 2500 deg 2 ) field. Although ambituous, a 
program like this is in the range of the near-future instru- 
mentation (e.g., VLT Survey Telescope, Visible & Infrared 
Survey Telescope for Astronomy) which encourages to per- 
form a comprehensive search for untriggered GRBs in the 
near future. In addition, spectroscopic follow-up observa- 
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tions of candidate orphans would be important to distin- 
guish between afterglows and other optical transients. 

The non-detection of an afterglow also provides a limit 
on the rate of other explosive events with similar fading 
behaviour. This includes events with minor or quashed 
high energy emission like failed GRBs (Huang et al. 2002) 
or so-called "dirty fireballs" (Dermer etal. 1999). 

6. Conclusion 

We presented the data reduction, analysis and results of an 
i?-band survey dedicated to slowly variable optical tran- 
sients. The survey strategy was designed specifically to 
search for afterglows of untriggered gamma-ray bursts. 
12 deg 2 were monitored in up to 25 nights down to a limit- 
ing magnitude of R=23. Throughout the survey, four pre- 
viously unknown transients were discovered. Based on the 
limited photometric data only putative classifications of 
the candidates could be obtained so far. The observations 
of three of the transients suggest them to be a cataclysmic 
variable, a flare star and a dwarf nova, respectively. 

The fourth transient appeared in a single image as a 
bright source slightly offset with respect to an underly- 
ing quiesent object. The decay of the source could not 
be followed due to the lack of observations, thus the or- 
gin of the transient is unresolved. However, the spatial 
association with the persistent source suggests an extra- 
galatic origin. The steep brightness increase of ~1.5 mag in 
~2 days appears a-typical for a supernova and thus makes 
the detected flaring source the best candidate for an or- 
phan afterglow in our survey. 

Simulations of the transient detection efficiency of the 
survey showed that the effective sky coverage was not 
sufficient to obtain a strong constraint on the collima- 
tion of the optically emitting GRB outflow from the non- 
detection of suitable counterparts. Nevertheless, we found 
that a similar programs like the one described in this pa- 
per are feasible to be performed in the near future. Limits 
on the collimation ratio of the X-ray to gamma-ray emit- 
ting regions were already obtained in the past and strong 
ratios (>8) were ruled out (Grindlay 1999; Greiner etal. 
2000). Tighter constraints on the optical emitting region 
will become available soon as well with the use of large 
dedicated surveys. 
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